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Abstract 



TV-dark-dark solitons in the generally coupled integrable NLS equations are de- 
rived by the KP-hierarchy reduction method. These solitons exist when nonlinear- 
ities are all defocusing, or both focusing and defocusing nonlinearities are mixed. 
When these solitons collide with each other, energies in both components of the 
solitons completely transmit through. This behavior contrasts collisions of bright- 
bright solitons in similar systems, where polarization rotation and soliton reflection 
can take place. It is also shown that in the mixed-nonlinearity case, two dark-dark 
solitons can form a stationary bound state. 

Keywords: Coupled nonlinear Schrodinger equations, KP hierarchy, dark-dark 
solitons, r function. 



1 Introduction 



In studies of nonlinear wave dynamics in physical systems, nonlinear Schrodinger (NLS)- 
type equations play a prominent role. It is known that a weakly nonlinear one-dimensional 
wave packet in a generic physical system is governed by the NLS equation pQ. Hence this 
equation appears frequently in nonlinear optics and water waves [2j [3j H] . Recently, it has 
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been shown that the nonlinear interaction of atoms in Bose-Einstein condensates is gov- 
erned by a NLS-type equation as well (called Gross- Pit aevskii equation in the literature) 
[5]. In these physical systems, the nonlinearity can be focusing or defocusing (i.e., the non- 
linear coefficient can be positive or negative), depending on the physical situations [I] or 
the types of atoms in Bose-Einstein condensates [5]. When two wave packets in a physical 
system or two types of atoms in Bose-Einstein condensates interact with each other, their 
interaction then is governed by two coupled NLS equations [2J [31 El El El El HOI ED] • The 
single NLS equation is exactly integrable [12J. It admits bright solitons in the focusing 
case, and dark solitons in the defocusing case. Its bright iV-soliton solutions were given 
in [12], and its dark iV-soliton solutions can be found in [13]. The coupled NLS equations 
are also integrable when the nonlinear coefficients have the same magnitudes [HI [151 EE] • 
In these integrable cases, if all nonlinear terms are of focusing type (i.e., the nonlinear 
coefficients are all positive), the coupled NLS equations are the focusing Manakov model 
which admits bright-bright solitons [14J. If all nonlinear terms are of defocusing type (i.e., 
the nonlinear coefficients are all negative), the coupled NLS equations are the defocusing 
Manakov model which admits bright-dark and dark-dark solitons [T71 [TBI 119] . If the focus- 
ing and defocusing nonlinearities are mixed (i.e., the nonlinear coefficients have opposite 
signs), these coupled NLS equations admit bright-bright solitons [16,120] and bright-dark 
solitons [21]. Existence of dark-dark solitons in this mixed case has not been investigated 
yet. 

Soliton interaction in these integrable generally coupled NLS equations is a fascinating 
subject. In the focusing Manakov model, an interesting phenomenon is that bright solitons 
change their polarizations (i.e. relative energy distributions among the two components) 
after collision [13] . In the coupled NLS equations with mixed nonlinearities, energy can 
also transfer from one soliton to another after collision [20]. In addition, solitons can be 
reflected off by each other as well [TO]. In the defocusing Manakov model, two bright-dark 
solitons can form a stationary bound state, a phenomenon which does not occur for scalar 
bright or dark solitons [17J. All these interesting interaction behaviors can be described 
by multi-soliton solutions in the underlying integrable system. In the focusing Manakov 
model, iV-bright-bright solitons were derived in [TJ] by the inverse scattering transform 
method. In the mixed-nonlinearity model, two- and three-bright-bright solitons and two- 
bright-dark solitons were derived in [201 EI] by the Hirota method, and iV-bright-bright 
solitons were derived in [16J by the Riemann-Hilbert method. In the defocusing Manakov 
model, iV-bright-dark solitons were derived in [17] , and degenerate two-dark-dark solitons 
were derived in [18] , both by the Hirota method. 

So far, progress on dark-dark solitons in the integrable generally coupled NLS equations 
is very limited. While dark-dark solitons in the defocusing Manakov model were derived 
in [18], we will show that their two- and higher-dark-dark solitons are actually degenerate 
and reducible to scalar dark solitons. In [19], the inverse scattering transform method was 
developed for dark solitons in the defocusing Manakov model. But, as we will show in this 
paper, their analysis can not yield general dark-dark solitons either due to their choices 
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of the boundary conditions. To date, general multi-dark-dark solitons in the coupled NLS 
equations have never been reported yet (to our knowledge). As we will see, these general 
multi-dark-dark solitons are not easy to obtain due to non-trivial parameter constraints 
which must be met. 

In this paper, we comprehensively analyze dark-dark solitons and their dynamics in 
the generally coupled integrable NLS equations. First, we show that these coupled NLS 
equations can be obtained as a reduction of the Kadomtsev-Petviashvili (KP) hierarchy. 
Then using T-function solutions of the KP hierarchy, we derive the general iV-dark-dark 
solitons in terms of Gram determinants. These dark-dark solitons exist in both the de- 
focusing Manakov model and the mixed-nonlinearity model. Recalling that bright-bright 
solitons exist in the mixed-nonlinearity model as well [TU [TS], we see that the coupled 
NLS equations with mixed nonlinearities are the rare integrable systems which admit 
both bright-bright and dark-dark solitons. The dark-dark solitons obtained previously 
in [T8J [19] for the defocusing Manakov model are only degenerate cases of our general 
solutions. Next, we analyze properties of these soliton solutions. For single dark-dark 
solitons, we show that the degrees of "darkness" in their two components are different 
in general. When two dark-dark solitons collide with each other, we show that energies 
in the two components of each soliton completely transmit through. This contrasts colli- 
sions of bright-bright solitons in these same equations, where polarization rotation, power 
transfer and soliton reflection can occur [TH [T6l [20] . Thus dark-dark solitons are much 
more robust than bright-bright solitons with regard to collision. In the case of mixed 
focusing and defocusing nonlinearities, an interesting phenomenon is that two dark-dark 
solitons can form a stationary bound state. This is the first report of dark-dark-soliton 
bound states in integrable systems. However, three or more dark-dark solitons can not 
form bound states, as we will show in this paper. 

We should mention that this KP-hierarchy reduction for deriving soliton solutions in 
integrable systems was first developed by the Kyoto school in the 1970s [22J. So far, 
this method has been applied to derive bright solitons in many equations such as NLS, 
modified KdV, Davey-Stewartson equations [231 EH E5j. This method has also been 
applied to derive iV-dark solitons in the defocusing NLS equation [25J. But this reduction 
for dark-dark solitons in the generally coupled NLS equations is more subtle and has 
never been done before. In this paper, we will derive general iV-dark-dark solitons by this 
KP-hierarchy reduction and the grace of deep use of determinant expressions. Compared 
to the inverse scattering transform method [19] and the Hirota method [18], our treatment 
is much more clean, and the solution formulae much more elegant and general. Thus, the 
KP-reduction method has a distinct advantage in derivations of dark-soliton solutions. 



3 



2 The TV- dark- dark solitons 



The generally coupled integrable NLS equations we investigate in this paper are 

iu t = u xx + (5\u\ 2 + e\v\ 2 )u, , 
iv t = v xx + (5\u\ 2 + e\v\ 2 )v, 

where 5 and e are real coefficients. This system is integrable [T4J [151 EES]- Through u and 
v scalings, the nonlinear coefficients 5 and e can be normalized to be ±1 without loss of 
generality When e = 5 = 1, this system is the focusing Manakov model which supports 
bright-bright solitons [14] . When e = 5 = —1, this system is the defocusing Manakov 
model which supports bright-dark and dark-dark solitons [UJ HHJ [19]. When e and 5 have 
opposite signs, the system exhibits mixed focusing and defocusing nonlinearities. In this 
case, these equations support bright-bright solitons [HHED], bright-dark solitons [2TJ, and 
dark-dark solitons (as we will see below). 

In this section, we derive the general formulae for iV-dark-dark solitons in the integrable 
coupled NLS system ([T]). The basic idea is to treat Eq. (jTJ) as a reduction of the KP 
hierarchy. Then dark solitons in Eq. ([T]) can be obtained from solutions of the KP 
hierarchy under this reduction. For this purpose, let us first review Gram-type solutions 
for equations in the KP hierarchy [26l [271 EH] • 



Lemma 1 Consider the following equations in the KP hierarchy 



{\D x D r - l)r(k) ■ r{k) = -r{k + l)r{k - 1) 
(D 2 X -D y + 2aD x )r(k + 1) • r(k) = 0, 



(2) 



where D is the Hirota derivative defined by 



d_ d_ 

dy dy' 



f(x,y)g(x',y')\x=x>, 



y=y > 



(3) 



a is a complex constant, k is an integer, and r(k) is a function of three independent 
variables (x,y,r). The Gram determinant solution r(k) of the above equations is given by 



r(k) =det 1 < i j< N (m ij (k) \ = m^k) 
where the matrix element rriij(k) satisfies 



l<i,j<N 



d x m i5 {k) = (pi(k)^j(k), 

d y mij(k) = (d x cpi(k))ipj(k) - cpi(k)(d x ipj(k)), 

d r mij{k) = -<fi(k - l)4jj{k + 1), 

rriij(k + 1) = mij(k) + (pi(k)ipj(k + 1), 



(4) 
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andifi(k) and ipj{k) are arbitrary functions satisfying 



dyipi(k) = dlipi(k), 
(fi(k + 1) = (d x - a)(fi(k), 
dyipj(k) = -d^j{k), 
tP J (k-l) = -(d x + a)ij J (k). 



Before proving this lemma, several remarks are in order. The first equation in ([2]) is 
the bilinear equation for the two-dimensional Toda lattice (see p. 984 of [29] and p. 4130 
of [2D]), and the second equation in ([2]) is the lowest-degree bilinear equation in the 
1st modified KP hierarchy (see p. 996 of |29j). Since the two-dimensional Toda lattice 
hierarchy and modified KP hierarchies are closely related to the (single-component) KP 
hierarchy, all these hierarchies will be called the KP hierarchy in this paper. Regarding 
the parameter a in the second equation in (J2J), it corresponds to the wave- number shift ko 
in [29] [see Eq. (10.3) there]. The bilinear equation with this parameter was not explicitly 
written down in [29J, but can be found in [30J [see Eq. (N-3) there]. This parameter can 
be formally removed by the Galilean transformation for y in ([2]). But for our purpose, it 
proves to be important to keep this parameter, as it will pave the way for the introduction 
of another similar parameter b in Lemma 2 later. In that case, a and b can not be removed 
simultaneously by the Galilean transformation, and they are essential for the construction 
of non-degenerate dark-dark solitons in the generally coupled NLS system ([T]). 

Proof of Lemma 1. By using (J4]) and f$5\j, we can verify that the derivatives and shifts 
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of the r function are expressed by the bordered determinants as follows 



d x r(k) 

dim 

d y r(k) 
d r r(k) 



+ 



rriijik) ifi(k) 

-4>j(k) 

m ij (k) d x ipi(k) 

-^(k) 

rriijik) d x ipi(k) 

-^(k) 

rriij{k) (pi(k-l) 

if>j{k + l) 



(d x d r - l)r(k) 



r{k + l) 

T{k - 1] - 
(d x + a)r(k + 1) = 



rriijik) <fi(k) 
-d x ^jik) 

rriijik) (fiiik) 
-d x ^jik) 



rriijik) Viik-l) <fi(k) 
ipjik + I) -1 
-ipjik) -1 

rriijik) <pi(k) 
-ipjik + 1) 1 



id x + afrik + 1 
id y + a 2 )rik + 1 



rriijik) ifi{k-l) 
1 

rriijik) d x (fi{k) 
—ipjik + 1) a 

rriijik) 8%<Pi{k) 



-ipjik + 1) a 2 

rriijik) d%(pi(k) 
-ipjik + I) a 2 



Here the bordered determinants are defined as 



rriij (fi 



mu "7,12 

m N1 rn N2 
-ipl ~ip2 



rriijik) d x piik) (p^k) 

-ipjik + 1) a 1 

-ipjik) 

rriijik) d x <fi(k) 

-^•(fc + l) a 1 

-ipjik) 



m 1JV </?i 



-^JV 



and so on. By using the Jacobi formula of determinants [26], we obtain the bilinear equa- 
tions ([2]) from the above expressions. □ 



Using Lemma 1, we can obtain solutions to a larger class of equations in the KP 
hierarchy below. 
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Lemma 2 Consider the following equations in the KP hierarchy, 



{\D x D r - l)r(k, I) ■ r(k, I) = -r(k + 1, l)r(k -1,1), 
(D 2 X -D y + 2aD x )r(k + 1,1)- r(k, I) = 0, 
{\D X D S - l)r(k, I) ■ r(k, I) = -r(k, I + l)r(k, I -I), 
(Dl -D y + 2bD x )r(k, l + l). r(k, I) = 0, 



(6) 



where a,b are complex constants, k,l are integers, and r(k,l) is a function of four inde- 
pendent variables (x,y,r,s). The solution r(k,l) to these equations is given by the Gram 
determinant 

r(k,l) = det 1 <i !j < N (rn ij (k,iyj = m^k,!) _ ^ (7) 
where the matrix element rriij(k, I) is defined by 



l<i,j<N 



mijik, I) = dj + ^-(fi(k, l)ipj(k, I), 
cpi(k, I) = {pi - a) k (pi - b) l e^, 



i \i 



ir 



qj +b< 



with 



& = PiX + v\y + ^z-j + ^z- b s + &o, 
Vj = QjX - q)y + + ^s + Vj0 , 

and Cij, Pi, qj, £j 0; r/j are complex constants. 



(9) 



It is noted that the system (JS)) is an expansion of the previous system (J2J) by adding a 
new pair of independent variables (s, /) to the previous pair (r, k). 

Proof. It is easy to see that functions m^{k, I), ifi(k, I) and i/jj(k, I) satisfy the following 
differential and difference rules, 

d x m,ij(k,l) — <fi(k, l)ipj{k, I), 

d y mij{k,l) = (d x ipi(k,l))ipj(k,l) - <p i (k,l)(d x '4) j (k,l)), 
d r mij(k, I) = —fi{k - 1, l)i>j(k + 1, /), 

mij{k + 1,1) = rriij(k, I) + <fi(k, l)^j(k + 1, 1), , . 

dytpi(k,l) = dg(pi(k,l), 1 

Lpi(k + 1,1) = {d x - a)(pi(k, I), 

d y ipj{k,l) = -d 2 x ijj(k,l), 

tpj(k -l,l) = -(d x + a)ipj(k,l). 

Then from Lemma 1, we can verify the first two bilinear equations in fl6]). The other 
two equations in can be obtained directly by replacing a, k, r as b, I, s in Eq. ([2]) of 
Lemma 1. □ 

Next, we perform a reduction to the bilinear system ([6]) in the KP hierarchy. Solutions 
to the reduced bilinear equations are given below. 
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Theorem 1 Assume that f is a real function of real x and t, and g, h are complex func- 
tions of real x and t, then the following bilinear equations 



(D 2 X + %| 2 + e\u\ 2 )f ■ f = 5\i2\ 2 gg + e\u\ 2 hh, 
(W t + D 2 X + 2icD x )g ■ f = 0, 
(iD t + D 2 X + 2idD x )h ■ f = 0, 



(11) 



where 5, e, c and d are real constants, /i and v are complex constants, and the overbar ' 
represents complex conjugate, admit the following solutions, 



f = 


Sij 


+ 


1 






Pl+Pj 


i 


9 = 


Sij 


+ 


1 


i Pt- 


ic 


Pl+Pj 


Pj+ic 


h = 




+ 


1 


1 Pt~ 


id 


Pi+pj 


V pj +id 



(12) 



where 

ij = /'/•'• + 'i>ji + bo, 

Pj are complex constants satisfying the constraint 



5\ji\ 



\Pj — ic\ 



+ 



€\V\ 



\pj — id\ 



-2, 



(13) 



(14) 



and £jo ar & arbitrary complex constants. 



Proof. In Lemma 2, if one assumes x, r, s are real, y, a, b are pure imaginary, k, I are 
integers, and qj = pj, r)j = ^j , Cji = c^, then we have 



(15) 



r]j = ij, m ji (k,l)=m ij (-k,-l), r(k, I) — r(—k, -I). 
Therefore, defining 

c l3 =5 lv Refe)>0, / = r(0,0), <? = t(1,0), h = r(0,l), (16) 

where 5ij is 1 when i — j and otherwise, then 

1 



/ = 



(0,0) 



Pi + Pj 



g = r(-l,0), h = r(0,-l), (17) 



and 



{\D x D r - l)f ■ f = -gg, 
{\D X D S - l)f ■ f = -hh, 
(D 2 X -D y + 2aD x )g ■ f = 0, 
(D 2 X -D y + 2bD x )h ■ f = 0. 



(18) 
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Under the above reduction, the solution fl7|) for r can be rewritten as 



r(k,l) 



1 / Pi~ a ^k, Pi-b- 



Pi+Pj Pj + a Pj + b' 



<\ i< 6 ^ + 



1 / Pi~ a ^k, Pt-b^i 



Pi+Pj Pj + a Pj + b' 



with 



& + 6 = (Pi + Pi)x + (p 2 - P 2 )y + (:— + —^-)r + (- - + ir — -) 5 + 6o 



>i - a Pi + a ft - & Pi + 6 



Thus if pi satisfies the constraint 



%| 2 (^— + ) + e M 2 (^T + -^t) = "2(p< + Pi), 

Pi - a Pi + a Pi - o Pi + 



i.e., 



+ 



6 1/ 



2. 



(Pi - a) + a) - 6) + b) 

then from Eqs. (JT5]) - (j2D]) . one gets 

(%| 2 <9 r + e|z/| 2 <9 s )r(A;, I) = -2d x r{k, I). 

Using / = t(0, 0), this equation gives 

S\^f r + e \u\ 2 f s = -2f x . 

Differentiation of ( 125|) with respect to x gives 

8\v\ 2 fxr + tW\ 2 fxs = ~Vxx- 

The first two equations of (1T81) are just 

fxrf fxfr / -99, 

fxsf - fxfs - f = -hh. 

So from Eqs. (1251) - (}2"oT) . we have 

Vxxf - VI + (%l 2 + e|H 2 )/ 2 = %l 2 ss + eW\ 2 hh, 

which is just 



Finally, denoting 



(D x + %| 2 + e\v\ 2 )f ■ f = 5\^gg + e\v\ 2 hh. 



y = it, a = ic, b = id, 



with t, c and d real, the second and third equations in (ITTj) and ( TT2"]) are obtained directly 
from Lemma 2, and the constraint (1141) is obtained directly from Eq. ( 12 ip . Theorem 1 is 
then proved. □ 

Now we transform the bilinear equations (TTTj) in Theorem 1 into a nonlinear form. To 

do so, we set 

9 - h 
u = fij, v = Uj, (30) 



where /, g, h satisfy Eq. (1111) . From ([301 . we have 

(D t g ■ f)/f = Zt/p, (D t h ■ f)/f = v t /v, 

(D x g ■ f)/f = u x /fi, (D x h ■ f)/f = v x /v, 
(D 2 x g ■ f)/f 2 = u xx /n + (u/»)(D 2 J ■ f)/f, (31) 
(D 2 x h ■ f)/f = v xx /u + (vM(Dlf ■ f)/f. 
The first bilinear equation in (II ip is 

D 2 J.f = -(%| 2 + e|H 2 )/ 2 + 5\rfgg + e\u\ 2 hh 

which can be further rewritten as 

(D 2 J ■ f)/f = -(%| 2 + e|H 2 ) + 5\u\ 2 + e\v\ 2 , (32) 

The second bilinear equation in (ITT]) is just 

{Dl + iD t + 2icD x )g-f = Q (3g) 

Substituting fl3Tl) into fl33|) . we have 

iu t + u xx + u{D 2 J ■ f)/f + 2icu x = 0. (34) 
In the same way, from the third bilinear equation in fTTTj) we have 

ivt + v xx + v{D 2 J ■ f)/f + 2idv x = 0. (35) 

Substituting ([32]) into dSD and ([35]), we get 

iu t + 2icu x + u xx + u[-5\n\ 2 - e\u\ 2 + 5\u\ 2 + e\v\ 2 } = 0, . . 

iv t + 2idv x + v xx + v{-5\/j>\ 2 - e\v\ 2 + 5\u\ 2 + e\v\ 2 } = 0. ^ ' 

Letting 

U = ue i[i-^\ 2 -M 2 +c 2 )t-cx]^ 
ft = ve i{(S\»\ 2 -e\v\ 2 +d' 2 )t-dx] 
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Eqs. f)36p are then transformed into 



iut + u xx + (S\u\ 2 + e\v \ 2 )u = 0, 
iv t + v xx + {5\u\ 2 + e\v\ 2 )v = 0, 



which has N-dark-dark soliton solutions as 

e i[cx+(S\, 

J[dx+(8\fi\ 2 +e\is\ 2 -d 2 )t] "hit 



U 
V 



.. A[cx+{5\^\ 2 +e\v\ 2 -c 2 )t\ gN_ 



ve 



(37) 



(3f 



/a 



with fN,9N,h N given by (fT2|) . Finally, taking t — > —t, Eqs. ( 1571) become the gener- 
ally coupled NLS equations (JTJ). Hence we immediately have the following theorem for 
solutions of Eq. (pQ). 



Theorem 2 The N-dark-dark soliton solutions for the generally coupled NLS equations 
([I]) are 



where 



n . . . r i\cx—(8\fj,\ 2 +e\u\ 2 —c 2 )t] G N 

a — fj,e Fn , 

i[dx-{&\^\ 2 +e\v\ 2 -d 2 )t\ Hn. 

Fn ' 



1 



V = 




F N = 


Sij + 


Gn = 


6ij - 


Hn = 


- 



(39) 



Pi+Pj 



NxN, 
J Pj-j£ e di+dj 



Pi+Pj Pj+lC 



1_ Pi-id c 0i+8j 



NxN 
NxN' 



(40) 



c,d are real constants, /i,u,pj,9jo are complex constants, and these constants satisfy the 
following constraints 



e\v\ 



\Pj 



ic\ 



\Pj ~id\ 2 



j = l,2,---,N. 



(41) 



These solitons are dark-dark solitons, i.e., both u and v components are dark solitons, 
because it is easy to verify that 



u _). ^ e i[cx-(S\ii\ 2 +€\i/\ 2 -c 2 )t+(f>±] 



v — > ve 



i[dx-(5\iJL\ 2 +e\v\ 2 -d 2 )t+x±} 



X — > ±00, 



(42) 



where <f>± and x± are phase constants. Thus the u and f solutions approach constant 
amplitudes and \v\ at large distances. When <5 > and e > 0, which correspond to self- 
focusing nonlinearities for both u and v components in Eqs. (0Q), the constraints (141 p can 
not be satisfied, thus dark-dark solitons can not exist as expected. When 5 < and e < 0, 
which correspond to self-defocusing nonlinearities for both u and v components, dark-dark 
solitons can exist as Ref. [19] shows. A new phenomenon revealed by Theorem 2 is that, 
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when S and e have opposite signs, which correspond to mixed focusing and defocusing 
nonlinearities in the u and v equations, the constraints (1411) can still be satisfied, hence 
dark-dark solitons can still exist. This phenomenon will be demonstrated in more detail 
in the next section. Interestingly, when 5 and e have opposite signs, Eqs. ([I]) also admit 
bright-bright solitons [16]. Thus Eqs. ([1]) with opposite signs of 5 and e are the rare 
equations which support both dark-dark and bright-bright solitons. 

The parameter constraints (j4"Tl) can be solved explicitly, so that solutions fl39|) can be 
expressed in terms of free parameters only. Let us write 



Pi 



cij +ibj, 



where aj and bj are the real and imaginary parts of pj. Then Eq. (I4ip becomes 



+ 



e\u\ 



a] + (bj - c) 2 a 2 + (bj - d) 



-2. 



Solving this equation, we find that a 2 can be obtained explicitly as 




-e\v\ 
2 1 1 



, _ C )2 _ {b . _ d) 2 + _ m 



1 , 

~e\v\ 

2 1 1 



5e|/z| 2 |z/| 2 



(43) 



(44) 



Here 5, e, /i, u, c, d and bj are all free parameters as long as the quantity under the square 
root of (I44p as well as the whole right hand side of (jUJ) are non-negative. If aj < 0, 
we will see that the soliton solution f[3"9j) would be singular. Thus in this paper, we will 
always take Oj > to avoid this singularity. 

We would like to make four remarks here. The first remark is on the above derivation 
of dark solitons through KP-hierarchy reduction. This derivation is non-trivial. To better 
understand it, we can split it into two parts. One part is the reduction of the bilinear 
equations (TTTT) of the generally coupled NLS equations (jTJ) from the KP-hierarchy equa- 
tions ([6]). The other part is the reduction of the soliton solutions to the bilinear equations 
(TTTjl from the r-solutions ((Zj) of the KP-hierarchy equations ([6]). In the first part, when 
we impose on the r- functions the conjugation constraint [see (TT5]) ] 



T(M)=T(-fc,-Z), 



and the linear constraint [see ( 122]) ] 



(5\fx\ 2 d r + e\u\ 2 d s )r(k, I) = -2d x r{k, I), 



(45) 



(46) 



and set 



/ = r(0,0), g = r(l,0), h = r(0,l), y = it, 



ic, 



id. 
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with t, c, d being real, then one can readily verify that the KP-hierarchy equations 
reduce to the bilinear equations (ITTj) of the coupled NLS equations ([T]). In the second 
part, in order for the r-functions (J7J to satisfy the conjugation constraint (I45j) . it is 
sufficient to require [see ffl~5]) ] 

rriji(k, I) = mij(-k, -/). (47) 
A sufficient condition for fj47|) to hold is that 

Cij = S i:j , qj = pj, i]j = rtjo = CjO) ( 48 ) 

x,r,s are real, and y,a,b are pure imaginary. These conditions are the same ones we 
imposed at the beginning of the proof of Theorem 1. Under these conditions, the r- 
solutions ([7]) of the KP-hierarchy equations fl6]) then reduce to the solutions fTT2|) for the 
bilinear equations fTTTj) of the coupled NLS equations ([T]). In order for the r-functions ([7]) 
to satisfy the linear constraint (146 p . by rewriting these r-functions as (1191) and inserting 
them into this linear constraint, we then get the parameter constraint f[2"Tj) . which is 
equivalent to the parameter constraint (l4Ti) in Theorem 1. This splitting of the earlier 
derivation of dark solitons into these two parts helps to clarify this derivation and make 
it more understandable. 

The second remark is on the solution form (j2HD of dark solitons in the generally coupled 
NLS equations ([T]). It is known that the NLS equation of focusing type is a reduction of 
the two-component KP hierarchy (see |29j, page 966 and 999), and the NLS equation of 
defocusing type is a reduction of the single-component KP hierarchy [25] . It is also known 
that solutions to the single-component KP hierarchy can be expressed as single Wronskians 
[2"B"| I3~T| I3"2"] , and solutions to the two-component KP hierarchy can be expressed as double 
Wronskians [33J. Thus iV-bright solitons in the focusing NLS equation can be expressed 
as double Wronskians |3H |35], and iV-dark solitons in the defocusing NLS equation can 
be expressed as single Wronskians [25J. These Wronskian solutions can also be expressed 
as Gram- type determinants [261 123 EE [36], [37J . For the vector generalization (JI)) of the 
NLS equation, in order to obtain its iV-bright-soliton solutions, one should increase the 
number of components, and take (CEJ) as a reduction of the three-component KP hierarchy. 
Thus iV-bright solitons in ([1]) can be expressed as three-component Wronskians (or the 
corresponding Gram- type determinants [26]). But to obtain iV-dark solitons in Eqs. 
([1]), one should increase copies of independent variables to (r, k) and {s,l) in the single- 
component KP hierarchy [see Eqs. (jBJ], thus iV-dark solitons in Eqs. (|T|) can still be 
expressed as single Wronskian (or the corresponding Gram determinant) as we have done 
above. 

The third remark we make is on comparison of the KP-hierarchy reduction method and 
the inverse scattering method for deriving dark-soliton solutions. As is well known, the 
inverse scattering method is another way to derive soliton solutions. For bright solitons, 
the inverse scattering method (or its modern Riemann-Hilbert formulation) is a powerful 
way to derive such solutions (see [381 [39] for instance). Recently, bright-bright iV-solitons 
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in a very general class of integrable coupled NLS equations were easily derived by this 
method [W\, and Eqs. (QQ) are special cases of such general equations. But for dark 
solitons, the inverse scattering method is more difficult due to non-vanishing boundary 
conditions, which create branch cuts and other related intricacies in the scattering process 
[13]. In [19], the inverse scattering transform analysis was developed for the defocusing 
Manakov equations [5 = e = — 1 in ([T])] with non- vanishing boundary conditions. But in 
their analysis, the boundary conditions 042p were taken such that c = d [see their equation 
(2.3)] (actually c = d — was taken there, but the case of c = d ^ can be reduced 
to the case of c = d = through Galilean transformation). When c = d, one can see 
from our general formula (1391) that u and v are simply proportional to each other, thus 
their inverse scattering analysis could only obtain degenerate dark-dark solitons which 
are reducible to scalar dark solitons in the defocusing NLS equation. In order to derive 
the more general dark-dark solitons (13 9 j) with c ^ d, the inverse scattering method would 
be even more complicated than that in [19]. Comparatively, the KP-hierarchy reduction 
method we used above is free of these difficulties, and is thus a simpler method for deriving 
dark-soliton solutions. 

Our last remark is on dark solitons in an even more general coupled NLS equations 

iu t = u xx + (5\u\ 2 + e\v\ 2 + juv + juv) u, , . 

ivt = v xx + (5|m| 2 + e\v | 2 + juv + juv) v , 

where S, e are real constants as in (QQ), and 7 is a complex constant. If 7 = 0, (I49p reduces 
to (P)). This more general coupled NLS system (j4"9~]l is also integrable. Its Lax pair as 
well as iV-bright-bright solitons are given in [16] . To explore dark-dark solitons in this 
system, we look for solutions with the following large- distance asymptotics [as in Q39]) ] 



where //, v are non- zero complex constants, and c, d, lo, k are real constants. Inserting this 
asymptotic solution into (}4"9~]) . we see that due to the 7-terms, Eqs. ( 14"9"|) can hold only 
if c = d, and u — k. Based on the previous solutions fl39|) . this would imply that the u 
and v components of dark-dark solitons in the general system (I49p must be proportional 
to each other, thus are equivalent to scalar dark solitons in the defocusing NLS equation. 
Except these trivial dark-dark solitons, Eqs. (I49p do not admit other dark-dark solitons 
of the form ( |50|) when 7 7^ 0. This is a dramatic difference between the cases of 7 = 
and 7 7^ in Eqs. (H9|) . Whether the general system (T49|) admits dark-dark solitons with 
background asymptotics different from (150]) is still unclear. 



3 Dynamics of dark solitons 

In what follows, we investigate the dynamics of single-dark-soliton and two-dark-soliton 
solutions in the generally coupled NLS equations (JT]). In the analysis of these solutions, 




X —7- — OO 



(50) 
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5 and e will be treated as arbitrary parameters. In the illustrations of solutions in the 
figures, we will pick 

6=1, e = -l, (51) 

which correspond to mixed focusing and defocusing nonlinearities. The reason for this 
choice is that dark solitons under such mixed nonlinearities have never been studied 
before. We will show that under these mixed nonlinearities, some novel phenomena (such 
as existence of two-dark-soliton bound states) would arise. Soliton dynamics under other 

6 and e values, such as in the defocusing Manakov equations where 6 = e = — 1, would 
also be briefly discussed when appropriate. 



3.1 Single dark solitons 

In order to get single dark solitons in Eqs. (CQ), we set iV = 1 in the formula (|39i) . After 
simple algebra, these single dark solitons can be written as 



u 



1 llp i[cx-(5\tJ,\ 2 +e\u\ 2 -c 2 )t} 



1 1 , p i[dx-{5\^i\ 2 +e\u\ 2 -d 2 )t] 
2 



i + yi + (yi - 1) tanh( gi+ ^ i+/?i ; 



l + zi + (zi- l)tanh(^p 



pi 



(52) 



(53) 



where 



9 X = p x x - ip(t + 0io, e Pl = l/(pi+pi), 
y 1 = («c-pi)/(ic + pi), z x = (id- px)/(id + pi), 
and fj,, u,pi,9io are complex constants satisfying 



6\n\ 



e\v\ 



\Pi ~ %c \ 



\P\ - id\ 2 



(54) 



or equivalently, a% is given by formula (1441) . where p\ = a% + ib\. This soliton would 
be singular if p\ + p\ < 0, i.e., a\ < 0. Thus we will require a\ > below to avoid 
singular solutions. It is easy to see that the intensity functions \u\ and \v\ of these dark 
solitons move at velocity — 2b\. In addition, they approach constant amplitudes and 
\v\ respectively as x — > ±oo. As x varies from — oo to +oo, the phases of the u and v 
components acquire shifts in the amount of 2(f>i and 2xi, where 



2-1 



(55) 



i.e., 20! and 2\\ are the phases of constants yi and Z\ respectively. Without loss of 
generality, we restrict — n < 2<pi, 2xi < 7r, i.e., — n/2 < <f>i,Xi — 7r /2- At the center of the 
soliton where 6i + 9% + p\ = 0, intensities of the two components are 



\U\ center = | A* I COS 01, \V\ center = \U\ COS^l- 



(56) 
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These center intensities are lower than the background intensities \fi\ and thus these 
solitons are dark solitons. Notice that the center intensities of the u and v solutions are 
controlled by their respective phase shifts 2<pi and 2xi, thus these phase shifts dictate 
how "dark" the center is. This general single dark-dark soliton fl52|) - fl53|) has been derived 
for the defocusing Manakov model before by the Hirota method in [TBI ITT] . In particular, 
a parameter constraint similar to (|54l) was given in [TB]. If c = d, then yi — zi, hence 
4>i = Xi- I n this case, the u and v components are proportional to each other, and have 
the same degrees of darkness at the center. This soliton is equivalent to a scalar dark 
soliton in the defocusing NLS equation, thus is degenerate. It is noted that the single- 
dark-dark soliton derived in [19] [see Eq. (5.8) there] corresponds to this degenerate type 
of dark-dark solitons. To illustrate, we take 

/i = 1, v = 2, , c = d = 0, Pi — vL5, 6*io = 0, (57) 

which satisfy the constraint (15"4"|) . Intensities of the solution (l52 |) - (|53p are displayed in Fig. 
UJa). This soliton is stationary, and both its u and v components are black (with zero 
intensity) at the soliton center. 

Non-degenerate single-dark-dark-solitons in Eqs. ([T]), however, are such that c ^ d. 
The u and v components in these solitons are not proportional to each other, thus are 
not reducible to scalar single dark solitons in the defocusing NLS equation. Since c ^ d, 
yi ^ zi, thus <pi Xi- This means that the u and v components in these non-degenerate 
solitons have different degrees of darkness at its center. To illustrate, we take 

^ = 1, u = 2, c = 0, d=0.5, pi = 1.0679, (58) 

which also satisfies the constraint (|54p . Here the pi value is obtained from the formula (|44j) 
with the plus sign and b\ = 0. Intensities of this soliton are displayed in Fig. QJJd). This 
soliton is also stationary. At its center, the u component is black, but the v component 
is only gray. This type of non-degenerate single dark-dark solitons in the coupled NLS 
system ([TJ has not been obtained before (to our knowledge). 

In the defocusing Manakov equations where 5 = e = — 1, their degenerate and non- 
degenerate single dark solitons qualitatively resemble those shown in Fig. [TJ and are thus 
not shown. 



3.2 Collision of two dark solitons 

Two-dark-soliton solutions in system (TjQ) correspond to iV = 2 in the general formula (1591 . 
In this case, we have 



u = lie 



F 2 (x,t) 



(59) 
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Figure 1: Single dark-dark solitons in Eqs. <Q with 5 = l,e = — 1: (a) a degenerate 
soliton with parameters ( 15"7"|) ; (b) a non-degenerate soliton with parameters (I5S|) . 

= i[rfx-(5|/,| 2 +e|H 2 -d 2 )t]^0M) ( 60 ) 

F 2 (x,t)' 1 ; 

where 

F 2 (x, t) = 1 + e 9l+e_1+pl + e e2+e_2+P2 + fe 9 i+^+^+^+pi+P2 ) ( 61 ) 

G 2 (x, t) = 1 + y ie ei+e " 1+Pl + y 2 e e2+e " 2+P2 + r yi y 2 e° l+§1+e2+ ^ +pl+p \ (62) 

F 2 (z, t) = 1 + Zl e ei+ ° 1+Pl + z 2 e e2+e " 2+p2 + rz lZ2 e ei+ ° 1+e2+ ° 2+pl+p2 , (63) 

0,- = Pj x - ip 2 jt + 9 j0 , e« = l/( Pj + Pj), (64) 

y . = (ic - Pj )/ (ic + pj-), ^ = (id - p 3 -)/ + Pi)> ( 65 ) 

r = 1 - (pi +pi)(p 2 +P2VIP1 +P2I 2 , (66) 

and /i, is,px,p2, 6*io, 6* 2 o are complex constants satisfying the constraint (1411 with j = 1, 2, 
or equivalently, cij is given by the formula (JHJ) , where pj = aj + ibj . 

In generic cases where Im(px) 7^ Im(p 2 ), these solutions describe the collision of two 
dark-dark solitons. To demonstrate these collisions, we take parameters 

/x = 1, v = 2, c = 0, d = 0.5, pi = 0.8426-0.2i, p 2 = 1.1801+0.2i, 1O = 9 20 = 0. (67) 

Here the real parts of pi and p 2 are obtained from the formula (jHJ) with the plus sign. 
The corresponding two dark-dark soliton solution (I59l) - (l60j) is shown in Fig. [2j We can 
see that after collision, the two dark solitons pass through each other without any change 
of shape and velocity in either of its two components. Hence the degrees of darkness in 
each soliton do not change after collision, which means that there is no energy transfer 
from one component to the other inside each soliton after collision. In addition, there 
is no energy transfer from one soliton to the other after collision either. This complete 
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transmission of dark solitons' energy in both its two components after collision occurs not 
only for 5 = 1 and e = — 1 as in Fig. [2J but also for all other S and e values. Thus it 
is a common phenomenon of the generally coupled NLS system ([1]). For instance, it also 
happens in the defocusing Manakov equations where 5 = e = — 1. 

This complete transmission of dark-dark solitons' energy in both its two components 
is a remarkable phenomenon, because it is in stark contrast with collisions of bright- 
bright solitons in the same coupled NLS system ([1]). Indeed, for bright-bright solitons in 
the focusing Manakov system (with 5 = e = 1), polarization rotations take place after 
collision, hence energy has transferred from one component to the other in each soliton 
[HJ . For bright-bright solitons in the more general coupled NLS system (such as 5 = 1 
and e = — 1 above), energy can also transfer from one soliton to another after collision 
[TB] . Thus collisions between bright-bright solitons and between dark-dark solitons in the 
coupled NLS system (JTJ are distinctly different. 

The reason for this complete energy transmission in all components in dark-soliton 
collisions is that the intensity profile of each dark-dark soliton is completely characterized 
by the background parameters /i, u, c, d and the soliton parameter pj [see Eqs. ( 152]) - ( 153]) ]. 
These background parameters are the same for both colliding solitons, and clearly do not 
change before and after collision. The soliton parameter pj corresponds to the spectral 
discrete eigenvalue in the inverse scattering transform method, and is a constant of motion 
throughout collision. Consequently, the intensity profile of each dark-dark soliton (in both 
u and v components) can not change before and after collision. This property indicates 
that dark solitons are more robust than bright solitons with regard to collision. The 
positions of dark solitons do shift after collision though, as can be seen clearly in Fig. [2] 
This position shift is always toward the soliton's moving direction, which is the same as 
collisions of bright solitons in the NLS equation |12j . 



4 Dark-dark-soliton bound states 



In studies of dark solitons, multi-dark-soliton bound states is an interesting subject. In the 
defocusing NLS equation, two dark solitons repel each other, thus can not form a bound 
state [ID]. In the defocusing Manakov model, multi-bright-dark-soliton bound states were 
reported in [17J. Some of those bound states are stationary, while the others are not. So 
far, multi-dark-dark-soliton bound states have never been reported in integrable systems. 
In a non-integrable system, namely, the second-harmonic-generation (SHG) system, two- 
dark-dark-soliton bound states do exist, as was reported in [41J. In this system, single 
dark-dark solitons with non-monotonic tails exist. When two such dark-dark solitons 
weakly overlap with each other and interact, their non-monotonic tails create local minima 
in the effective interaction potential, hence the two dark-dark solitons can form stationary 
bound states. In addition, some of these bound states are stable jH]. 
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Figure 2: Collision of two dark-dark solitons in Eqs. ([I]) with 5 = 1, e = —1 and parame- 
ters (j67p . The upper row shows the (x, t) evolution, and the lower row shows the intensity 
profiles before and after collision: t — — 10 (solid); t = 10 (dashed). 

In this section, we show that in the generally coupled NLS system ([T]), when both S 
and e are negative, i.e., all nonlinear it ies are defocusing (the defocusing Manakov model), 
multi-dark-dark-soliton bound states can not exist. But for mixed focusing and defocusing 
nonlinearities, where 5 and e have opposite signs, two-dark-dark-soliton bound states 
do exist and are stationary. To our knowledge, this is the first report of multi-dark- 
dark-soliton bound states in integrable systems. Properties and physical origins of these 
stationary bound states in the mixed-nonlinearity model ([I]) are quite different from the 
stationary bright-dark-soliton bound states in the defocusing Manakov model [41] and 
stationary dark-dark-soliton bound states in the non-integrable SHG model |41j . as we 
will explain later in this section. 

To obtain dark-dark-soliton bound states, the two dark solitons in the solution flo^]) - 
(160]) should have the same velocity, i.e., Im(pi) = Im(p 2 ) (or b\ = 6 2 )? so that the two 
constituent dark solitons can stay together for all times. In order for this to happen, 
two different (positive) values ai and a 2 from Eq. (T43j) must exist for the same values of 
b\ = &2- When S and e are both negative, where the nonlinearities are all defocusing, this is 
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not possible. The reason is that when 5 < and e < 0, the function on the left side of Eq. 
(I43p is an increasing function of a?. Thus for this function to reach the value level of —2 
on the right side of Eq. (j43[) . there is at most one a| solution, hence at most one positive 
dj value. This means that when nonlinearities are all defocusing (i.e., the defocusing 
Manakov model with S = e = — 1), there are no multi-dark-dark-soliton bound states. 
However, when 5 and e have opposite signs, where focusing and defocusing nonlinearities 
are mixed, the function on the left side of Eq. ( l4"3"j) may become non-monotone in a?, 
hence it becomes possible for Eq. f l4"3"j) to admit two different positive values a\ and a 2 
for the same values of b\ = b 2 (see below). In the formula (|44|) . these different a\ and 
a 2 values correspond to the plus and minus signs respectively. In this case, two-dark- 
dark-soliton bound states would exist, and this is a new phenomenon in the coupled NLS 
equations (JT|) under mixed focusing and defocusing nonlinearities. Physically, these results 
on bound states in Eqs. ([TJ can be heuristically understood as follows. We know that in 
the scalar defocusing NLS equation, two dark solitons repel each other. In the coupled 
NLS system ([TJ, if 5 and e are both negative, all nonlinearities are defocusing, hence two 
dark-dark solitons still repel each other, and no bound states can be formed. However, if 
S and e have opposite signs, parts of the nonlinear terms are focusing, and the other parts 
defocusing. While the defocusing terms repel two dark solitons, the focusing terms do 
just the opposite, which is to attract two dark solitons. Thus, when these repulsive and 
attractive forces balance each other, two dark-dark solitons then can form a stationary 
bound state. This physical mechanism for the existence of dark-dark-soliton bound states 
is quite different from that in the SHG model jJTJ (see earlier text). 

Next we examine these two-dark-dark-soliton bound states in more detail. Through 
Galilean transformation (i.e., in the moving coordinate system with this common velocity), 
this common velocity can be reduced to zero. Hence pi and p 2 become real parameters. 
In this case, it is easy to see that this bound state becomes 



7/ _ . lf 4\cx-(S\u\*+eM 2 -<?)t] G 2 (x) 

v _ h , c i\dx~(8\u\ 2 +eM 2 -d 2 )t] H 2 {x) 

F 2 (x) ' 



(68) 



where 



pJ x ) = 1 + e 2 Pi x + 2a i+' 3 i _|_ g2p2^+2o2+P2 _|_ re 2pix+2p 2 x+2a 1 +2a 2 +pi+P2 ^gg^j 
QJx) = 1 + y ie 2 P^ x+2a ^+P^ _|_ y 26 2p2X+2a 2 +P2 _|_ r y^y^ e 2p 1 x+2p 2 x+2a 1 +2a 2 +pi+p2 
H 2 (x) = 1 + Zie ^+ 2a ^ + Z2e lP2X+2a 2 +p 2 + rz ^ e 2p lX +2p 2 x+2a 1 +2a 2 +p 1 +p 2 ^ ^ 

aj = Re(9jo), and Pj,Vj,Zj,r are as given in Eqs. (I64j) - (l66l) . Notice that functions 
F 2 , G 2 and H 2 are time- independent, thus this bound state is actually stationary. This is 
analogous to certain bright-dark-soliton bound states in the defocusing Manakov model 
[IT] and dark-dark-soliton bound states in the SHG model [41J. An important feature 
of these present bound states is that, as x moves from —00 to +00, these states acquire 
non-zero phase shifts. Indeed, it is easy to see from the above solution formula that the 
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phase shifts of the u and v components are 

w-phase shift = 20i + 2(p 2 , f-phase shift = 2xi + 2\2, (72) 

where 2<pj and 2\j are the phases of i/j and Zj respectively. In other words, the total 
phase shifts of the bound state are equal to the sum of the individual phase shifts of the 
two constituent dark solitons, which are non-zero in general. This contrasts stationary 
bright-dark-soliton bound states in the defocusing Manakov model [T7] and dark-dark- 
soliton bound states in the SHG model [41] . where phase shifts of the dark components 
across the soliton are all zero. 

To demonstrate these stationary two-dark-soliton bound states, we take parameters 

fi = 1, v = 2, c = 0, d = 0.5, pi = 1.0679, p 2 = 0.3311, ot x = a 2 = 0. (73) 

Here p± and p 2 are obtained from the formula (Pill) with b% = b 2 = 0. The corresponding 
bound state is displayed in Fig. [3] (upper row). In this bound state, the w-component is 
double-dipped (i.e., has a double hole), signifying this is a two-soliton bound state, while 
the -y-component is single-dipped. By adjusting a\ and a 2 values, we can obtain bound 
states where both u and v components are double-dipped. For instance, when we take 
Oil = — ot-2 = 2 instead of zero in (T73]) . we get such a bound state which is shown in the 
lower row of Fig. [31 For both bound states, the total phase shift of the w-component is 
zero, and the total phase shift of the f-component is 3.4355, as can be calculated from 
formula (T72]). 

From the above analytical formulae and Fig. [31 we can see that these stationary two- 
dark-soliton bound states have six free parameters, /i, u, c, d, a,\ and a 2 [the positive p\ 
and p 2 values are determined from formula (|44p by setting b\ = b 2 = 0]. The first 
four parameters characterize the background intensities and phase gradients, while the 
parameters ot\ and a 2 control the positions of the two dark solitons. 

The above dark-soliton bound states in the integrable coupled NLS system ([I]) possess 
properties which are very different from those in dark-soliton bound states in the non- 
integrable SHG model [4T]. First, the bound states in the SHG model are formed by 
identical dark solitons (see Fig. 4 in |4T]). but the bound states in the coupled NLS 
system are formed by different dark solitons since p± ^ p 2 (see lower row of Fig. [3] in 
this paper). Second, the bound states in the SHG model have zero phase shifts from 
one end to the other, but the phase shifts of bound states in the coupled NLS system 
([T]) are non-zero in general [see Eqs. (172"]) ]. Thirdly, the bound states in the SHG model 
have non-zero binding energy, hence can be stable against perturbations [41]. But the 
bound states in the present coupled NLS system have zero binding energy. Thus under 
perturbations, the two constituent dark solitons in these bound states generically will 
split apart, analogously to bright-soliton bound states in the focusing NLS equation. 

At this point, one may wonder if three- and higher-dark-dark-soliton bound states exist 
in the coupled NLS system (TjQ). It turns out that such bound states can not exist. The 
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Figure 3: Two examples of two-dark-soliton bound states in Eqs. ([I]) with S = 1, e = — 1. 
Upper row: bound state with parameters (1731) ; lower row: bound state with parameters 
(173")) except that a\ = —a 2 = 2. The left two panels show (x,t) evolution of \u\ and \v\ 
components, and the right panel shows the stationary intensity profiles. 



reason is that, in a bound state, velocities of all constituent solitons must be the same, 
i.e., all bj [i.e., Im(pj)] must be the same. In order for three- and higher-dark-soliton 
bound states to exist, formula ( 144"|) must give at least three distinct positive solutions aj 
for the same bj value. This is clearly impossible, since formula (jHJ) can give at most two 
distinct positive aj values when the plus and minus signs are taken. Consequently, three- 
and higher-dark-dark-soliton bound states can not exist in Eqs. (1T1) . Note that in the 
defocusing Manakov model, non-stationary three and higher bright-dark-soliton bound 
states exist, but stationary three and higher bright-dark-soliton bound states do not [TT] : 
while in the non-integrable SHG system, stationary three and higher dark-dark-soliton 
bound states do exist [41] . 



22 



5 Summary and discussion 



In this paper, we have investigated dark-dark solitons in the integrable generally coupled 
NLS system ([T]). By reducing the Gram- type solution of the KP hierarchy, we derived 
the general iV-dark-dark solitons in this system. We showed that the dark-dark solitons 
derived previously in the literature are only degenerate cases of these general soliton 
solutions. We have also shown that when these solitons collide with each other, energies 
in both components of the solitons completely transmit through. This behavior contrasts 
bright-bright solitons in this system, where polarization rotation and soliton reflection 
can occur after collision. In addition, we have shown that when focusing and defocusing 
nonlinearities are mixed, two dark-dark solitons can form a stationary bound state. These 
results will be useful for many physical subjects such as nonlinear optics, water waves and 
Bose-Einstein condensates, where the coupled NLS equations often arise. 

The dark-dark solitons obtained in this paper for the generally coupled NLS system ([T]) 
are useful for other purposes as well. For instance, it is known that from dark solitons 
of the defocusing NLS equation, one can obtain homoclinic solutions of the focusing 
NLS equation through simple variable transformations |42]. Thus, from these dark-dark 
solitons in this paper, we can obtain homoclinic solutions to these generally coupled NLS 
equations (pQ). Since solutions near homoclinic orbits often exhibit chaotic dynamics [12], 
the homoclinic solutions for the generally coupled NLS equations (jTJ) then can serve as 
the starting point to understand chaotic behaviors in these systems. 

Lastly, we would like to mention that iV-bright-bright and iV-bright-dark solitons in the 
coupled NLS equations ([T]) can also be obtained by the KP-hierarchy reduction method. 
But those reductions will be different from the ones in this paper for dark-dark solitons, 
and will be left for future studies. 
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